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Abstract. We analyse the anisotropy of the cosmic microwave background (CMB) 
for the Poincaré dodecahedron which is an example for a multi-connected spherical 
universe. We compare the temperature correlation function and the angular power 
spectrum for the Poincaré dodecahedral universe with the first-year WMAP data 
and find that this multi-connected universe can explain the surprisingly low CMB 
anisotropy on large scales found by WMAP provided that the total energy density 
parameter Qtot is in the range 1.016... 1.020. The ensemble average over the primordial 
perturbations is assumed to be the scale-invariant Harrison-Zel’dovich spectrum. The 
circles-in-the-sky signature is studied and it is found that the signal of the six pairs of 
matched circles could be missed by current analyses of CMB sky maps. 


PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 


1. Introduction 


The question of the global geometry of the Universe, i.e. of its spatial curvature, 
topology and thus its shape, is a very fundamental one. While the problem of the 
spatial curvature is discussed in all textbooks on cosmology, mainly in connection with 
the total mass/energy densities in the Universe and the various inflationary scenarios, 
the possibility of a multi-connected Universe is only rarely mentioned. This is surprising 
in view of the fact that it was already realized before Einstein’s seminal first paper on 
cosmology [I] that the three-dimensional space of astronomy might be not only non- 
Euclidean, but also multi-connected according to a given Clifford-Klein space form (see 
e.g. |). If the Universe is assumed to be simply connected, the cosmological principle 
implies that it is not only locally, but also globally isotropic and homogeneous. This 
means that all spatial points are geometrically equivalent. In the case of a multi- 
connected Universe, however, the cosmological principle holds in general only on the 
universal covering space of constant curvature, i.e. on E3, S? and H? for a flat, positively 
and negatively curved Universe, respectively. While the Einstein gravitational field 
equations still hold in a multi-connected Universe, the global structure of the Universe 
at large scales will be more complicated. 

The discovery of the temperature fluctuations ôT of the cosmic microwave 
background radiation (CMB) by COBE in 1992 [3] and the detailed measurements 
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by WMAP [4] and by other groups have led to a renewed interest in the question of the 
global geometry of the Universe (see [5] [6] for reviews). 

The most promising signatures to detect a possible multi-connected spatial 
structure of the Universe are the strange suppression of the CMB quadrupole and 
octopole, and of the temperature correlation function first observed by COBE [7] and 
nicely confirmed by WMAP [f], and the so-called circles-in-the-sky-signature proposed 
in [B]. 

Until about 1999, all observational data were consistent with the fact that the 
total energy density Etot (at the present epoch) amounted only to ca. 30% of the 


critical energy density Ecrit = 3Hoe ie. Qiot = Etot/Ecrit ~ 0.3 < 1, and thus 
strongly indicated that the geometry of the Universe is hyperbolic. Consequently, several 
groups performed detailed studies of compact and non-compact models of the Universe 
possessing hyperbolic geometry. Although later observations on the magnitude redshift 
relation of the supernovae of type Ia [9] [0] provided evidence for a large dark energy 
component, and the determination of the first acoustic peak in the CMB anisotropy 
by TOCO [fi] 12), BOOMERanG [13] and MAXIMA-1 [4 [5] indicated a nearly 
flat Universe, i.e. Qst œ~ 1, the available data are still compatible with the spatial 
geometry of the Universe being hyperbolic. (For a discussion of hyperbolic universes, 
see [0 LA MA 1 £0 20 22 eA eA 5) A 271.) 

The concordance model of cosmology assumes an exactly (spatially) flat Universe 
with the topology of £3 and, furthermore, that the dark energy is given by a positive 
cosmological constant A, i.e. Qa := AG = 1 — Qmat — Orad With Qmat = Opar + Qeam; 
where the various parameters denote the present values of the baryonic (bar), cold 
dark matter (cdm), matter (mat) and radiation (rad) energy densities in units of Ecrit 
(ACDM model). In the following, we shall refer to three variants of the concordance 
model presented by the WMAP team on the Legacy Archive for Microwave Background 
Data Analysis (LAMBDA) Web site |http://lambda.gstc. nasa.gov] see also [28]. These 
models give a good overall fit to the CMB anisotropy on small and medium scales, but 
there remains a strange discrepancy at large scales as first observed by COBE [7] and 
later substantiated by WMAP [4]. 

Let us expand the temperature fluctuations ÔT (À) of the microwave sky, where the 
dipole contribution has been subtracted, into real spherical harmonics Yim (ft) on S?, 


MOE SST aim MON (1) 


l=2 m=-l 
where ñ denotes the unit vector in the direction from which the photons arrive. (Note 
that the monopole and dipole terms, l = 0,1, are not included in the sum (J).) From 
the real expansion coefficients aj, one forms the multipole moments 


Ci = gah È Cm) 2) 


m=— 
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and the angular power spectrum 


(1+1 
ÔT? := vt Lo. (3) 
T 


The average (...) in (2) denotes an ensemble average over the primordial perturbations 
to be discussed in section M| respectively an ensemble average over the universal 
observers. The temperature two-point correlation function CV) is defined as C(V) := 
(6T (n)dT (A'Y with À- nr’ = cos V, which can be computed from the multipole moments 
(2) under the assumption of statistical isotropy as 


1 OO 
CV) ~ in (21+1)C,;P(cosv) . (4) 
l=2 
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Figure 1. C(V) from first year WMAP data (solid curve) in comparison to the three 


concordance models provided by the WMAP team at http://lambda.gsfc.nasa.go 


model pl assumes a power law fit to the WMAP, CBI and ACBAR data, model ri 
assumes a running index fit to the WMAP, CBI and ACBAR data, and finally, model 
bf assumes a running index fit to the WMAP, CBI and ACBAR data and additionally 
takes the large scale data from 2dF and Ly-alpha into account. 


In figure []we show C (9) as solid curve as measured by WMAP [4]. The important 
conclusion to be drawn from figure[JJis that the temperature correlation function displays 
very weak correlations at wide angles, 60° < J < 160°. The data are compared with 
three different versions of the concordance model obtained by the WMAP team as best 
fits of their data combined with various other data. The model pl (dotted curve) assumes 
a power law fit, i.e. with no running index, to the WMAP [4], CBI [29] 80) and ACBAR 
[BI] data, the model ri (dashed curve) assumes a running index fit to the WMAP, CBI 
and ACBAR data, and finally, the model bf (dashed-dotted curve) assumes again a 
running index fit to the WMAP, CBI and ACBAR data and additionally takes the 
large scale structure data from 2dF and Lyman-a [83] into account. It is seen from 
figure [] that the difference between the three concordance models is not considerable. 
One makes, however, the important observation that the three concordance models 
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using the best-fit values for the cosmological parameters as obtained by WMAP do not 
reproduce the experimentally observed suppression of power at wide angles, leaving 
us with “the mystery of the missing fluctuations” [84]. In the case of the largest 
angles, 160° < 3 < 180°, where the WMAP curve in figure [I] shows even negative 
temperature correlations, one might argue that the measurements in the nearly back- 
to-back configuration are very subtle due to the need to apply a Galactic cut to the 
data and, furthermore, due to geometrical reasons, since the expectation values of C(v) 
are computed from a much smaller number of pixels, and thus the discrepancy at these 
angles should not be overemphasized. There is, however, also a discrepancy between 
the three concordance models and the WMAP data at small angles, 9 < 40°, which is 
clearly visible in figure [I] 

Depending on certain priors, the WMAP team reported for the total energy 
density Qto = 1.02 + 0.02 together with Qpar = 0.044 + 0.004, Qmat = 0.27 + 0.04, 
and h = 0.71*0:03 for the present day reduced Hubble constant (the errors give the 1o- 
deviation uncertainties). Taken at their face value, these parameters hint to a positively 
curved Universe. Recently Luminet et al. [85] studied the Poincaré dodecahedral space 
(for details, see section 2) which is one of the well-known space forms with constant 
positive curvature (see also BE [37] [34]). The authors of ref. [B5] computed the CMB 
multipoles for l = 2,3 and 4, fitted the overall normalization factor to match the WMAP 
data at l = 4, and then examined the prediction for l = 2 and l = 3 as a function of 
Qiot- For eot = 1.013 they found a strong suppression of the power at l = 2 and weak 
suppression at l = 3 in agreement with the WMAP data. Since the eigenfunctions of the 
Poincaré dodecahedron are not known analytically, they must be computed numerically 
[38] [39]. In ref. [85], only the first three modes with wave number 6 = 13,21 and 
25 (comprising in total 59 eigenfunctions) have been used, which in turn restricted 
the discussion to the multipoles l < 4. There thus remains the question about how 
this extremely low wave number cut-off affects the prediction of the multipoles, since 
experience shows that increasing the cut-off usually enhances the integrated Sachs- Wolfe 
contribution. 

In this paper, we present a thorough discussion of the CMB anisotropy for the 
dodecahedral space topology which is based on the computation of the first 10521 
eigenfunctions corresponding to the large wave number cut-off 6 = 155. Taking within 
the tight-coupling approximation not only the ordinary, but also the integrated Sachs- 
Wolfe and also the Doppler effect into account, we are able to predict sufficiently 
many multipole moments such that a detailed comparison of the dodecahedral model 
with the WMAP data can be performed. We find that the temperature correlation 
function for the dodecahedral universe possesses very weak correlations at large scales 
in nice agreement with the WMAP data for Qot in the range 1.016...1.020. Besides 
the suppression of large-scale fluctuations, the dodecahedral model predicts 6 pairs of 
matching circles in the sky possessing an angular radius of 40°...50° for the slight 
curvature of space corresponding to the cited range of Qtot. 

We begin in section] by describing the spherical space forms and, in particular, the 
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Poincaré dodecahedron, and discuss the general properties of the vibrational modes. In 
section B] we outline our numerical method to determine the eigenmodes numerically. 
In section [4] we calculate the angular power spectrum and the temperature correlation 
function for the dodecahedral space and compare the predictions with the WMAP data. 
As a quantitative test, we study the S statistic introduced in [28]. For this statistic, it is 
found that nice agreement with the WMAP data is obtained for Qy,, = 1.016... 1.020. 
In section E] we study the circles-in-the-sky signature using a quantity © introduced in 
[8] (called S in [8]). In the last section, there is a discussion and summary. 


2. The spherical space forms and their vibrational modes 


The three-dimensional spaces M? of constant positive curvature K = +1 (spherical 
spaces) were classified by 1932 [40] 41] and are given by the quotient M? = S3/T of the 
3-sphere S? under the action of a discrete fixed-point free subgroup [ C SO(4) of the 
isometries of S3. All these manifolds are compact and are, apart from their universal 
covering S’, multi-connected. 

The unit 3-sphere S? is defined as the three-dimensional hypersurface of the 
four-dimensional unit ball in flat four-dimensional space with Cartesian coordinates 


(w,z,y,z) € Rf, ie. S is given by the set of points which satisfy the condition 


w?+a?+4+y?+ 22 = 1. The spatial distance |dz| between neighbouring points is 
given by 

dz? = dw? + dr? + dy? +d2? , (5) 
and the full four-dimensional space-time line element reads 

ds’ = êdt? — R°(t)dz? = a(n) [dn — dz’°] . (6) 


Here R(t) = a(n) denotes the cosmic scale factor as a function of cosmic time t 
respectively conformal time 7 (dn = cdt/R(t)). The metric (@) has the most general 
form assuming that the Universe is homogeneous and isotropic (cosmological principle). 

The condition w? +z? +y?+2z? = 1 can be used to eliminate one spatial coordinate, 
w say, and one is left with the three space variables (x,y,z). In the following, we shall 
also use the three angular variables (,6,) defined by 


w= cosy , x = sinxsin@cos¢d , 
y =sinysindsing , z = sinycosé , (7) 


withO<y<7,0<60<7,0<¢< 2r. It is seen from (Y) that (x,y,z) describe for 
a fixed value of y an S?-sphere with radius sin y and volume VS’ (x) = 1(2y — sin2y), 


where 0 and ¢ play the familiar role of polar and azimuth angle. The S?-sphere with 
radius sin y can be considered as a cross section of the unit S*-sphere. Varying x 
between 0 and 7, one obtains an infinite sequence of cross sections (=2-spheres) of S3, 
whose radius grows from 0 to 1 (for 0 < x < 4), and then shrinks again to zero (for 
= <x <7), which provides a useful way to visualize the 3-sphere. 
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With @ and (J, the line element on S? can be rewritten as 
dz? = dx? + sin? x (dé? + sin” 6 dd’) (8) 
which has to be inserted into eq. (6). 

For a study of the large-scale anisotropies in the CMB produced by scalar 
perturbations, we need the (regular) eigenmodes of the covariant Laplacian A on S’, 
i.e. the solutions of the Helmholtz equation (A+ ES*)YS? = 0 describing the vibrations 
on S®. Using the coordinates (y,6,¢) and the metric (8) and separating the variables, 
the normalized solutions read [42] [43] [44] 47] (8 € Njl = 0,1,...,8— l;m = 
—l,-1+1,...,0) 

ie. ~ 
Viim?) = Ret(x) Ym(0, $) . (9) 
Here Yim(0,@) are real spherical harmonics on the unit sphere S?, and the “radial 
functions” Rg(x) are given by 


Rai(x) = Ap (sin x) C5ti(cosx) , (10) 


where C(x) are the Gegenbauer polynomials (which can also be expressed in terms of 
the associated Legendre functions P?(x)). The normalization factor in (£0) 


1, / (8-1-1)! 

Ay = 2ta ly a 11 

k a (8+)! E 
follows from the orthonormality relations 

f du(E) Veim(@) Yrm (E) = Sap" Sw Sram’ (12) 
S3 

i dx sin? x Rai(x) Rau(x) = sp (13) 
0 


The volume element on S? reads (with the standard expression dQ = sin dôd on S?) 
du(Z) = sin? xdxdQ . (14) 
Notice that the radial functions Rg)(x) are even or odd (depending on whether 6 —/—1 


is an even or odd integer) about the equator y = 5, since they satisfy the periodic 


boundary condition 


Raia —x) = (-1)°"?! Rax) . (15) 
Since S? is compact, the vibrational modes are discrete, 
3 


where 6 € N plays the rôle of the wave number. The fundamental mode of S3 has an 
eigenvalue of ES’ = 0 corresponding to the constant wave function US) = 1/VV*, 
where VS = f 33 du(Z) = 27° is the volume of the unit sphere S?. Since the eigenvalues 
of S? do not depend on l and m, the modes are degenerate, and the multiplicity of the 
mode £ is given by 


rO = S\(214+1) = 8. (17) 
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If NS? (k) denotes the number of the vibrational modes on S3 with wave number smaller 
than or equal to k € N, we obtain 

NS (k) = nps M4 iR ER O(k? 

#) = Dr) = bak + ek = ak + OF) 

in agreement with Weyl’s asymptotic law for k — oo. 

We now turn to the multi-connected spherical spaces S3/T, T C SO(4), possessing 
the volume V(S3/T) = V°’/N, where N denotes the order of the group T. (For 
reviews, see KO [41] [48] [49] [50}). To define the discrete fixed-point free subgroups 
T C SO(4) of the isometries of S°, it is advantageous to work with Hamilton quaternions 
q := w + ri + yj + zij, (w, x,y,z) € Rt, which are spanned by the 4 basis quaternions 
{1,i,j,ij}, with the multiplication defined by i? = j? = —1, ij = —ji plus the property 
that i and j commute with every real number. The multiplication of quaternions is 
associative, but not commutative. If g* := w — xi — yj — zij denotes the quaternion 
conjugate to q, the norm of q is defined by |q|? := qq* = w? + 2? +y? + 2°. Quaternions 
with norm 1 are called unit quaternions. For a unit quaternion q, there exists an inverse 
q7} which is simply given by q7! = q*. It is obvious that the unit 3-sphere S? can 
be identified with the multiplicative group of unit quaternions. The distance d(q1, q2) 
between any two points qı and q2 on S? is given by cos d(q, q2) = W1Wet2122+Y1y2+2122 
and is a point-pair invariant, i.e. d(yq, yq2) = d(q, q2) Vy E€ SO(4). 

The group SO(4) is isomorphic to S? x S*/{-+(1, 1)}, the two factors corresponding 
to the left and right group actions. In the following, we are only interested in subgroups 
T € SO(4) which lead to homogeneous manifolds S?/T. Such groups possess only 
so-called right-handed Clifford translations. A quaternion y € I corresponds to a 
Clifford translation, if it translates all points q1,q2 € S? by the same distance y, 
i.e. d(qı, yq1) = d(q2,7¥GQ2) = x. The right-handed Clifford transformations yp € T 
in S? are defined by left-multiplication of an arbitrary unit quaternion q € S? by 
Yk = Ap + bki + Ckj + dyij, 1. e. 

q = q= (ak + bri + ckj + drij) - (w + zi + yj + zij) 
akw — bkt — Cky — dkz 


1 


= ( 
(bkw + akz — dky + ckz 
( 


) 
) 
CkW + dpx + apy — bkz) j (18) 
+ (dw — cpu + bky + akz) ij 
The right-handed Clifford translation (I8) acts as a right-handed corkscrew fixed-point 
free rotation of S°. 

For our numerical computations of the eigenfunctions, it is convenient to represent 
the action (£8) of a group element yp € T by the following orthogonal 4 x 4 matrix My 


Ak —bk — Ck —dk 
bk Qk —d;, Ck 
Mk = 19 


dk Cy by, ak 
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which acts from the left in an obvious way on the basis quaternions (1, i,j, ij). 
The following groups lead to homogeneous manifolds S?/T MO AD [48] [49]: 
e The cyclic groups Z,, of order n (n > 1). 
e The binary dihedral groups D*, of order 4m (m > 2). 
e The binary tetrahedral group T* of order 24. 
e The binary octahedral group O* of order 48. 
e The binary icosahedral group J* of order 120. 
Let us consider the binary icosahedral group J* in more detail consisting of N = 120 
group elements. They are generated by the two right-handed Clifford translations [50] 
Nig > yq = jq = —y + zi + a oy 
Yya:q => Vq 
with the unit quaternion %2 := $ + si + al corresponding to the angles (V2, 02, ¢2) = 
(era) 
domain) of the binary polyhedral space D := S?/I* is a regular dodecahedron and is 


arctang), where o = vör is the golden ratio. The fundamental cell (Dirichlet 


known as the Poincaré dodecahedral space, which is made of 12 pentagons, 30 edges and 
20 corners. 120 dodecahedra tessellate the unit 3-sphere S*. The shortest translation 


distance in D is given by d(q, %24) = X2 = f, while 7 translates all points already to 
the much larger distance xı = 5- 

The vibrations on the general homogeneous spherical 3-spaces MÌ = S3/T are 
determined again by the Helmholtz equation on S?, but now with the periodic boundary 
conditions corresponding to the defining fundamental cells (polyhedra). Thus the 
eigenfunctions U“"(#) on M? can be expanded into the eigenfunctions US’ (Z): 

6-1 ı 


wee) = E Y Em M?) VSB), (20) 


with EM? = 62-1, 8 € N, i = 1,...,r"(6). It is important to notice that for 
a given manifold M? the wave numbers 6 do not take all values in N. The allowed 
B values together with their multiplicities r (8) are, however, explicitly known BI]. 


The full information of the non-trivial topology of M? is contained in the real expansion 
coefficients {€%),,(M°)} which by virtue of (I2) are given by 


i > 3 > 34 > 
Shin ME) = | AME) Win) WHOM) (21) 
They satisfy the normalization condition 


B-1 Į 
YY Emm) = N (22) 


l=0 m=-l 
due to the orthonormality relation 
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Summation of the relation i over the degeneracy index i yields 


22 


m? (g 


3 (Ebim (M J’ = N r’ (B) . (24) 


r 


In section B] we shall determine the coefficients Ebim (M?) numerically up to some wave 
number cut-off Bmax. Here we state only the following Conjecture for homogeneous space 
forms (0 <1 < 6-1) 


I r™?’(B) 


M*(B 
wri de (S5im(M D? =v) (25) 


m=—l i=l 


for which a numerical check will be discussed below. Let us also note that in our 
application to the CMB anisotropy, we only need the modes with 6 > 3 (if they exist), 
since the values 3 = 1,2 correspond to modes which are pure gauge terms [52]. 


3. Numerical determination of the eigenmodes 


For the homogeneous 3-manifolds M? = S/T, which we are considering here, the 
eigenvalues of the Helmholtz equation are given by Be = B?-—1, B EN, and even 
their multiplicities r^? (B ) are explicitly known [BI]. As already mentioned in section B] 
the wave number ( does not assume all values in N, and thus the multiplicity r(B) 
refers only to the 8 values which appear in the spectrum. E. g. in the case of the Poincaré 
dodecahedral space D, there exist no eigenmodes for even (-values, and, furthermore, 
also among the odd values 6 € 2N + 1 there are finitely many gaps, i.e. @ values for 
which no eigenmode exists. Explicitly, one has [51] 


E? = p? -1 (26) 
with 6 € {1, 13, 21, 25, 31, 33, 37, 41, 43, 45, 49, 51, 53, 55,57} U {2n + 1,n > 30} and 


r= ale] E E  @ 


The fact that the eigenvalues respectively the wave numbers 6 for a given 
manifold M? are explicitly known, facilitates greatly the numerical computation of the 
eigenfunctions yig ), since it saves us from the time-consuming numerical search 
for the allowed -values as it has to be carried out, e.g. for hyperbolic manifolds in 
the case of a negatively curved universe [53], [54]. Furthermore, since the expansion (20) 
in terms of the eigenfunctions S? involves for a given wave number only a finite basis 
of dimension 8° (see eq. (C), the number of the real coefficients {€%),,(M°)} is also 
restricted to 8? for a given eigenfunction with wave number 8 and degeneracy index i. 

The eigenfunctions have to satisfy the fundamental periodicity condition 


i 34 
Tila) = VG , Yges , Vuer , (28) 


where 7, denotes a group element of the group T C SO(4) which defines the manifold 
M? = S3/T. We will use the condition (28) in a collocation algorithm. Thus, the 
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condition (28) will be imposed on L > £6? randomly distributed points q on S°. 
To generate such random points, we use instead of the coordinates (w,2,y,z) the 
coordinates (v, œa, y) defined by 

w = yl-— 2w cosa 
x = vV2v cosy 
y = yw siny 

z = vyl-2vsina 
with v € [0,5] and a, € [0,27]. The line element () takes then the form 

dv? 

2u(1 — 2v) 


(29) 


di? = + (1—2v) da? + 2udy* _ , 

which gives the volume element du = dv dady. The random points q on S? are now 
generated by randomly choosing the coordinates (v,a,y) and then calculating from 
(29) the coordinates (w, x,y,z) of the corresponding unit quaternion q. In this way, we 
generate L random points {q1, q,---,qz}. An arbitrary group element y € I maps the 
point qr to the point g := Ykqr on S?, where in the following computation the group 
elements yk run over ['\{1}. The action of yẹ will be represented by the matrix Mk, 
see eqs. [3 19). After having generated in this manner L mirror points {%, @,.--, dz}, 
one imposes the condition (28) at all those L points. Inserting the expansion (20), one 


finally obtains the following over-determined system of equations (r = 1,..., L > 6?) 
for the 8? coefficients {&)),,(M°)}: 
B-1 1 
Yo [Yin (Ge) — im) Eim M’) = 0 (30) 
1=0 m=-l 


The system (80) is for a given -value of the form i, By; = 0, r = 1,..., L, where 
the (L x 8?)-matrix B = B() is given by the square bracket in eq. (80) and is thus 
completely determined in terms of the eigenfunctions on the covering space S. The 
vector € = (8) has the 8? components é; = €g1m(M2), i.e. s denotes the multi-index 
{00,1 — 1,10,11,...,8— 18 — 1}. The system (80) is numerically solved by using the 
singular value decomposition. One then obtains rM* (8) linearly independent solutions 
for the expansion coefficients {€%),,(M°)} which are distinguished by the index 7. In a 
final step, the coefficients are normalized according to the condition (22). 

Using this method, we have computed the expansion coefficients for T = J*, i.e. 
for the Poincaré dodecahedral space D = S?/I*, for 8 < 155 comprising the first 10521 
eigenfunctions. The normalized coefficients have been used to check the conjecture (25) 
which has been found to hold within a numerical accuracy of 13 digits. In addition, 
we have computed the coefficients {&%),,(M°)} and checked the conjecture (25) for 
B < 33 for some cyclic groups, some dihedral groups, the binary tetrahedral and binary 
octahedral group and have found again excellent numerical agreement. However, for 
the inhomogeneous lens spaces L(12,5) and L(72,17), we have found that the 
conjecture (25) does not hold. We thus conclude that the conjecture (25) is only valid 
for the homogeneous 3-spaces. 
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4. The angular power spectrum ôT? and the correlation function C(v) 


In order to determine the angular power spectrum ôT? and the temperature correlation 
ô 


function C(v), one has to compute the temperature fluctuation (ñ). This is done 
using the Sachs-Wolfe formula which is given within the tight-coupling approximation 
by (c= 1) 


T ia) = for iGo + au -+ rane GORRO] 


no p- 
+2 fer | ay a 
NSLS n 


where the k-integration has to be replaced by a summation over ( in the case of a discrete 


N=NSLS 


Ve(r(),9,6) , (31) 


spectrum. (Here nsus denotes the conformal time at recombination corresponding to a 
redshift zsrs = 1089.) ue are the eigenfunctions on M3, and ®;(7) are the expansion 
coefficients of the metric perturbation ®, see eq. (82) below. (7(7),0,¢) denote the 
spherical coordinates of the photon path in the direction ñ, where we now assume that 
the observer is at the origin of the coordinate system, i.e. at (Xobs, Jobs; Pobs) = (0, 0, 0). 
Here, ôg (n) is the expansion coefficient of the relative perturbation in the radiation 
component, and Vz (n) is the expansion coefficient of the spatial covariant divergence 
of the velocity field of the tightly coupled radiation-baryon components. The method by 
which the quantities occurring in (I) are numerically computed, is described in detail 
in Section 2 of [26]. 

The ordinary Sachs-Wolfe (SW) contribution to the temperature fluctuation dT (À) 
is a combination of the gravitational redshift due to the gravitational potential 
®(n,7,0,¢) at the surface of last scattering (SLS) and the intrinsic temperature 
fluctuation tô, due to the imposed entropic initial conditions. Since the SW 
contribution, i.e. the first two terms in the brackets in (BI) (see also eq. (83) below), 
gives a dominant contribution to the CMB anisotropies at large scales, we first present 
an exact analytic expression for the expansion coefficients {am }, the multipole moments 
{Ci} and the temperature correlation function C(V) for the Poincaré dodecahedral 
universe taking only the scalar perturbations into account due to the SW effect. 

For an energy-momentum tensor Ty with Tj; = 0 for i # j and i,j = 1,2,3, the 
line element can in conformal Newtonian gauge be written as [52] 


ds? = a(n) [(1+2®)dn? — (1 — 28)jdz |°] . (32) 


Denoting by n = n(@,¢) the unit vector in the direction in which the photons are 
observed, the SW formula reads 


STEW (À 1 
TA = O(nsis, Tss, 9,6) + gor iss, Ters, 9, 9) (33) 
1 
= 3 (Mss, Tsis: 9, >) (34) 


with Tss := No — nsus. (Here no denotes the conformal time at the present epoch.) The 
SW contribution (83) can be approximated as in (84) for modes well outside the horizon 
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27_| The approximation 


at the time nss of last scattering, i. e. 
(84) is used in the following derivation of the analytic expressions for CSW, eq. (EI), and 
CSW (9), eq. (44). In the actual numerical computation, we use (BI) as discussed above. 
For the dodecahedron D, the metric perturbation can be written as an expansion in the 
eigenfunctions of D 


®(n,T, 0, d) => SE aa TE Oo) . (35) 
B>13 i=1 
Here the prime in the summation over the modes 8 indicates that 8 € 2N + 1 does 
not assume all values. The functions ©%(7) determine the time-evolution and will be 
factorized 


®5() = ©5(0) gln) (36) 
with gg(0) = 1. The functions gg(7) do not depend on i, since the associated differential 
equation depends only on the eigenvalue ER which is independent of i. The initial values 
',(0) play the rôle of the primordial fluctuation amplitudes and are Gaussian random 
variables with zero expectation value and covariance 


(D500) ®5(0)) = dae dir Po(3) (37) 


Here Pə(6) denotes the primordial power spectrum which determines the weight by 
which the primordial modes £ are excited, on average. The average (...) in (87) denotes 
an ensemble average over the primordial perturbations. In the following, we shall assume 
the scale-invariant Harrison-Zel’dovich spectrum 
a 
P9(8) = A) 

The only free parameter in (88) is the normalization constant œ which will be determined 
from the CMB data. 

Inserting the expansion (85) into the Sachs-Wolfe formula (84) and using the 
eigenfunctions of the dodecahedron (see eqs. (20), @) and (ij), we obtain for the 
expansion coefficients of the SW contribution (l > 2) 


(38) 


r? (8) 
= XO 4 (sis) Ebim (D) Rai(tsrs) (39) 
oa i=1 


which in turn yields with the definition (2) and the Gaussian hypothesis (B7) the mean 
values of the multipole moments (l > 2) 


l 


j? 
->> Ps(8) 93 (nts) yo pp a (Esim (P Ra(Tsıs) - (40) 


ET m=—l i=1 
At this point we can use our Conjecture with N = 120 and arrive at the final form 
40 «~ r?(8)Ps(8 
cwo) = SI HORO gisis) Ral) (41) 


B>l 
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The last expression shows that the lowest multipoles are suppressed due to the discrete 
spectrum of the vibrational modes of D. Since the lowest modes are at 6 = 13, 21, 25 
and 31, the mode sum starts at 6 = 13 for the multipoles l = 2,...,12, while for the 
multipoles l = 13,...,20 it starts only at @ = 21, for the multipoles l = 21,...,24 only 
at 6 = 25, and for the multipoles l = 25,...,30 only at 6 = 31. 

From (4) we can compute the SW contribution to the correlation function (4) for 
the Poincaré dodecahedron D 


GY@) = i g3(nsis) C(B, v) ; (42) 


p213 
where the angular dependence is completely given by (P;(cos?) = g” *(cos V)) 


8—1 
C(8,9) = X (21 + 1) Rå, (Tss) Pi(cos ¥) 


= (21 + 1) A3 (sin Tsis)” [C31 (C08 Tars)] Ci? (cos v) (43) 

[=2 
after having inserted the radial function (IQ). With the help of the addition theorem 
for the Gegenbauer polynomials (see p. 223 in [55]), one obtains for the mean value of 
the SW contribution to the correlation function for the Poincaré dodecahedron universe 


CSW (9) = de r?(8)Pe(8) 92 (nsts) sin(37(V)) 


7 s 
3r = B sin 7(V) 
1 A 3 
= m E o cosù , (44) 
where the relation C3_, (cosy) = ien has been used. The angle y = (V) € [0, 27815] 
is given by 
(0) = arccos ( cos? Tests + sin? Tsg cos V) , (45) 


and the monopole and the dipole moment read 


asw _ 80 yr! r?(6)Pe(b) sin( rsrs) | 
Co = T2 — p i sls) es) (46) 
pe 
Gc = — sin?(Tsts) X A Blas) [C3_2(cos Tsis)” . (47) 
6213 


In figure P] we show as a dashed curve the evaluation of eq. (44) for Qot = 1.020. 
The solid curve shows the full result including the integrated Sachs-Wolfe and Doppler 
contribution. The grey band corresponds to the lo deviation computed from 500 
simulations. Although the curve corresponding to the ordinary Sachs-Wolfe contribution 
lies within the lo band of the full result, the other two contributions cannot be neglected 
as discussed in more detail below. 

In our computations of the multipole moments C; for a given primordial realization, 
we use all eigenmodes with @ < 155. Above this cut-off, we use in the mode summation 
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Figure 2. The SW contribution (44) of the temperature correlation function C(v) is 
shown as a dashed curve for the dodecahedral space D for Qtot = 1.020. The solid curve 
shows the full result including the integrated Sachs-Wolfe and Doppler contribution. 
The lo deviation is shown as a grey band. 


up to 8 = 1501 the conjecture (25) which describes the mean asymptotic behaviour. If 
only the modes below 6 = 155 are used, the multipoles can only be calculated up to 
l ~ 15. It is thus important for the calculation of C(V) that the conjecture allows us to 


include higher multipoles. 

Let us now come to the discussion of the CMB anisotropy in the Poincaré 
dodecahedral space, where, for simplicity, we assume that the dark energy is given 
by a cosmological constant with Qa := hrot — Qraa — Qmat. We have calculated the CMB 
anisotropy for different values of the cosmological parameters. Here we present only the 
0.28 and h = 0.70. 


results for Qpar = 0.046, Qmat 


10000 | 
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Figure 3. The values of the S(p) statistic (48) for the median values of C (9) for the 
dodecahedral space for p = 60° (full curve) and p = 20° (dotted curve) in dependence 
on Mot (h = 0.70, Amat = 0.28 and Qpar = 0.046). 


To quantify the observed surprisingly low CMB anisotropy [4] on large scales, the 
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S statistic 
cos p 
S(p) -J |C(0)|? dcos 0 (48) 


Ì 
is discussed for the first year WMAP data in [28] for p = 60°, and it is found that only 
0.3% of the simulations based on the concordance model ri have lower values of S(60°) 
than the observed value S(60°) = 1644. Somewhat higher values of S(60°) are obtained 
using other statistical methods and other sky masks in [56], but they are nevertheless 
surprisingly low. In figure]we show the values of $(60°) (full curve) and S(20°) (dotted 
curve) using the median values of C (V) for the dodecahedral space D in dependence on 
Qitot- One observes that the models with Qiot = 1.017...1.020 give the lowest values for 
the large scale anisotropy. 

Here and in the following calculations, the overall amplitude of the CMB anisotropy, 
i.e. the normalization constant a in (88), is fitted to the C; values of the first year WMAP 
data in the range l € [20, 45]. (Thus, the values cannot be directly compared to the ones 
in [85], where the amplitudes are scaled such that they match the C4 value of WMAP 
exactly.) 
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Figure 4. The temperature correlation function C (V) for the dodecahedral topology 
for four values of Qtot for h = 0.70, Qmat = 0.28 and Qpar = 0.046, the WMAP first 
year observation and the pl-concordance model (see description in the text). 


In figure{]the temperature correlation function Cv) is shown for four values of Qot- 


The WMAP curve obtained from the LAMBDA home page|http://lambda. gsfc.nasa.go 
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is shown as a full curve. The mean value of C(v) for the pl-concordance model from 
the same home page is displayed as dotted curve. The lo deviation computed from 
3000 HEALPix simulations is shown as a shaded region. Finally, the mean value of 
C (V) for the dodecahedral topology is shown as a dashed curved and the corresponding 
lo deviation is the grey region obtained from 500 simulations. In figure 4h), the value 
preferred in [B5], i.e. Qot = 1.013, is presented. One observes a general low anisotropy 
on large scales, however, a more than lo deviation from the WMAP values occurs in 
the range V € [60°, 90°]. The WMAP data are much better described by models with 
Qtot = 1.018 and ot = 1.019 (figures Mp) and ØẸ)) as is also revealed by the better 
values for the S statistic (compare figure B). Only on the largest scales J > 170°, the 
models in figures Ap) and ff) do not reproduce the large negative correlation obtained 
from the WMAP data. But without this blemish, the dodecahedral space describes in 
the considered Qiot range the observations down to the smallest scales, since the WMAP 
curve lies completely within the grey lo region. Finally, figure ØH) shows a model with 
Qtot = 1.026, where the S statistic already deteriorates and consequently the agreement 
with the WMAP curve is also worse. 
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Figure 5. The angular power spectrum ôT? for the dodecahedral topology (circles) 
for four values of Otot for h = 0.70, Qmat = 0.28 and Opa, = 0.046. The WMAP data 
are shown as diamonds with error bars, and the pl-concordance model is represented 
by the full curve. 


In figures Bh) to BH), the angular power spectrum ôT? = I(l + 1)C;/(27) is shown 
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for the same values of Qot as in figure] The WMAP data are represented by diamonds 
together with the lo error bars which do not include the cosmic variance. The mean 
values of ôT? of the pl-concordance model (full curve) do not show the suppression of 
power at low values of l, but in contrast, increase there. The mean values of ôT? for the 
dodecahedral space exhibit the suppression of power as observed by WMAP. A good 
agreement is observed exactly for those values of hot favored by the S statistic and 
the correlation function C (V). This behaviour is also revealed by figure [6] which shows 
the dependence on Miot of the mean values of ôT, ôT? and ôT? for the dodecahedral 
space D. The lowest values of the quadrupole ôT? are obtained near Q44 = 1.017. 
Furthermore, at Qot = 1.017 the next multipole ôT? matches the value obtained by 
WMAP almost perfectly. 


1.02 
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Figure 6. The mean values of ôT? (solid curve), ôT? (dashed curve) and ôT? (dotted 
curve) are shown for the dodecahedral space D as a function of Qtot. The lower 
horizontal line indicates the 6T? value of WMAP and the upper one the corresponding 
value of 6T?. 


The probability distribution of ôT? and ôT? obtained from sky simulations is for 
the pl-concordance model much broader than for the dodecahedral space. To emphasize 
this fact, figure [] displays the distribution of ôT? and ôT? for Qo = 1.018. The 
corresponding WMAP values are indicated by the vertical lines together with the lo 
error (vertical dark grey band). The shaded histogram shows the wide distribution 
of 3000 HEALPix simulations for the pl-concordance model. In contrast, a much more 
pronounced peak is revealed by the corresponding distribution of the dodecahedral space 
(light grey histogram). Furthermore, the peaks for ôT? and ôT? lie both within the 10 
error of the WMAP observation. 


5. The circles-in-the-sky signature 


A multi-connected space can reveal itself by the so-called circles-in-the-sky signature 
[8]. Consider an observer located at a given point 7 who measures the fluctuations of 
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Figure 7. The probability distribution of 57? and 6T? in [uK?] for Qtot = 1.018, 
h = 0.70, Omat = 0.28 and Qpar = 0.046 for the dodecahedral space (light grey 
histogram). The wide distribution of the pl-concordance model is shown as the shaded 
histogram. The WMAP values together with their lo error are indicated by the vertical 
line and the vertical grey band, respectively. 


the CMB. The contribution to the CMB fluctuations due to the ordinary Sachs-Wolfe 
(SW) effect arises on a sphere around this observer having the radius Tss, i.e. at the 
so-called surface of last scattering (SLS) w.r.t. this observer. In the universal cover 
there are copies of this observer which have their own SLS’s which are centered at the 
copies of the observer. Those copies whose distance to the observer at 7 are smaller 
than 27,3 possess SLS’s which intersect the SLS of the observer sitting at 7. Since the 
intersection of two spheres is a circle, both the observer and its copy observe along this 
circle the same temperature fluctuations due to the ordinary Sachs-Wolfe effect. Since 
the observer and its copy have to be identified, there are two matching circles having 
the same radius on the sky and having an identical ordinary Sachs-Wolfe contribution. 
The integrated Sachs-Wolfe (ISW) contribution arises on the photon path from the SLS 
to the observer, see eq. BIJ. Since the two paths are not identified, their contribution is 
different on the matched circles. The two Doppler contributions are also different, since 
their velocity projections are different. 

In the case of a flat toroidal universe, these different contributions to the 
temperature anisotropy are investigated in [58] with respect to the circles-in-the-sky 
signature, and it is found that the Doppler and the integrated Sachs- Wolfe contributions 
degrade the signature. That these two contributions play an important role also in the 
case of the dodecahedral topology, is demonstrated in figure [8] It shows the three 
competing contributions as computed from eq. (BI) for the simply connected S? (figure 
Bh)) and for the dodecahedral space D (figure Bb)). Here we use the conjecture (25) 
up to 6 = 4001. For the SW contribution, see eq. (EI). The most obvious difference 
between D and S? is the suppression of power for small multipoles ôT, in the case of 
D in comparison with S3. In the case of the dodecahedral space D, the loss of power 
is so severe that for the quadrupole, the ISW dominates over the ordinary Sachs- Wolfe 
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contribution. A further distinction between S? and D is that ôT, oscillates as a function 
of l in the case of the multi-connected space even for the mean values of ôT, which 
are considered here. With the exception of the smallest multipoles, the ordinary Sachs- 
Wolfe contribution dominates up to l X 30. However, even in this range the Doppler and 
the ISW contribution provide a significant fraction of the total temperature fluctuation. 
The ISW is most important for l < 10, whereas the Doppler contribution increases for 
l = 10 to even larger values than the ISW, such that both contributions give together a 
large perturbation to the pure circles-in-the-sky signature. 
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Figure 8. The three contributions to ôT; are shown for the simply connected S? (left) 
and the dodecahedral space D (right) calculated in the tight-coupling approximation 
for Mot = 1.019, h = 0.70, Amat = 0.28 and Qa, = 0.046. 


To emphasize this point, figure D] shows the temperature fluctuation along two 
matched circles parameterized by the angle ¢ defined in (49) below. One observes 
in figure Dh) that the total temperature fluctuation dT is not matched perfectly well. 
However, some rough similarities between both ôT curves are nevertheless visible. 
Whether such similarities occur sufficiently frequently by chance in the simply connected 
S? model, such that the circles-in-the-sky signature is swamped, is discussed below. The 
different contributions to ôT are separately shown in figures 9p), DF) and DH). Even the 
ordinary Sachs-Wolfe contribution (figure [9p)) does not provide a perfect match which 
is due to the fact that the dipole contribution has been subtracted from the sky map 
analogously as it is done on the observational side. This destroys the perfect agreement 
due to the SW contribution. Conversely, if one would know the topological structure of 
the Universe, this would offer an opportunity to determine the dipole contribution due 
to the primordial fluctuations which is usually superseded by the Doppler shift due to 
our local motion. The Doppler and the ISW contributions are completely different for 
the two circles as seen in figures DF) and OH), respectively, for reasons described above. 

As just discussed, for cosmological models near Qot œ~ 1.02, the ordinary Sachs- 
Wolfe contribution dominates for 1 < 30. The value l = 30 corresponds to a scale 
0 ~ 1802 = 6°. Thus the circles are blurred on scales below ~ 6°. On larger scales, the 
integrated Sachs-Wolfe and the Doppler contribution lead to a modulation of matched 
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Figure 9. The temperature fluctuation ôT is shown along two matched circles for 
Otot = 1.019, h = 0.70, Qmat = 0.28 and Qpar = 0.046. Panel a) displays ôT of the 
total temperature fluctuation along the two circles, obtained from the ordinary Sachs- 
Wolfe (SW), Doppler (Dop) and integrated Sachs-Wolfe (ISW) contribution, which are 
shown separately in panels b), c) and d), respectively. 


circle structure. 
For a quantitative search of matched circles in microwave sky maps, the quantity 


5p) = CITED O + 0) 
(TECE) + 6T5(¢)) 

is introduced in [B]. (In B], © is called S, a name we avoid in the following in order 
to not confuse the reader with the previously discussed © statistic.) Here ôT,(¢) and 


(49) 


ôT (Q) are the temperature fluctuations along two circles on the SLS with the same 
radius, and () := + - dọ. The angle p describes the relative phase between the 
temperature fluctuations on the two circles and the plus/minus sign in the nominator 
allows for orientable as well as for non-orientable manifolds. The © variable takes 
the value X = +1 in the case of a perfect correlation and © = —1 for anticorrelated 
temperature fluctuations. 

In [59] the © test is applied to the first year data from WMAP and it is found 
that Umax = max,{h(p)} does not show any spikes which would reveal a matched circle 
pair. The search is carried out for nearly back to back circles having a separation larger 
than 170° and a radius larger than 25°. For 1.0091 < Mso < 1.0252 (for h = 0.70 
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Figure 10. The radius of the six matched circles is plotted as a function of OQtot. 


and Qmat = 0.28) the Poincaré dodecahedron has in total 6 pairs of circles which are 
back to back, i.e. are separated by 180°. The radius of these circles increases with 
increasing ot as shown in figure [0] Above Qiot = 1.0252 there are at first 16 pairs, a 
number which increases further with increasing 2)i.4, and below Qo, = 1.0091 there are 
none. The conclusion in is thus that the topology of the Poincaré dodecahedron is 
excluded. However, in a hint towards six pairs of circles having a radius of 11 + 1° 
is found which would correspond to Qi. = 1.0094 for our choice of h and Qynat. 

Now we would like to discuss whether the Doppler and ISW contributions are 
significantly large enough in order to contaminate the SW contribution such that the 
circles-in-the-sky signature is hidden. We generate 500 sky maps for the dodecahedral 
space D using all modes up to 6 < 155 and compute for each of the six pairs of matched 
circles the value of Xmax. In this way we get 3000 values for Xmax belonging to circle pairs 
including the SW, Doppler and ISW contribution. In addition, we obtain 3000 values 
for 4(p) where p is not varied but instead is fixed to the value p = 36° determined by 
the topology of D. Due to the Doppler and the ISW contribution, the values of Xmax are 
usually larger than the correct values X(p) for the matched pair of circles having p = 36°. 
This emphasizes that the contamination by the Doppler and the ISW contribution is so 
large that the Xmax statistic does in general not find the correct value of p. Furthermore, 
we compute also the corresponding values of Xmax for 500 simulations for the simply 
connected S? universe with the same cut-off in 8 and with the ring positions found in 
the dodecahedral space D. Of course, the Xmax values of S? do not correspond to a 
matched circle pair since S? has none. 

In figure [I| we compare the probability distribution of Xmax for S? and D for 
the two values Qio = 1.011 and Qi., = 1.020. These values are chosen in order to 
demonstrate the effect due to the size of the radii of the circles which are 25.4° and 
49.8° for OQio¢ = 1.011 and rot = 1.020, respectively. The probability distributions 
shift to smaller values of Xmax with increasing radii. In figure [2] we again show the 
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probability distribution of Xmax for S? as in figure [I] but display for D the probability 
distribution of U(p) with p = 36°. It is seen that the distribution for D has been 
shifted to smaller values of £, and thus the overlap of the distributions for D and S° 
has increased. 


4+ Qot = 1.011 J 4 
PO max) PO max) 


0 G h pe í Pa, 0 HS a BA i 

0.0 0.2 0.4 0.6 0.8 Xmax 1.0 0.0 0.2 0.4 0.6 0.8 Snax 1.0 
Figure 11. The probability distribution of Umax for S? (shaded histogram) and 
D (grey histogram) for Qtot = 1.011 in panel a) having a radius of 25.4° and for 
Otot = 1.020 in panel b) having a radius of 49.8°. The other cosmological parameters 
are h = 0.70, Qmat = 0.28 and Opar = 0.046 as usual. 
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Figure 12. The probability distribution of Xmax for S? (shaded histogram) as in 
figure [I] and of £(p) for D with p = 36° as determined by the topology of D (grey 
histogram) for Qtot = 1.011 in panel a) having a radius of 25.4° and for Otot = 1.020 in 
panel b) having a radius of 49.8°. The other cosmological parameters are as in figure 


m 


Analysing a sky map which possesses no pairs of circles, one would obtain for 
randomly chosen circle pairs a distribution corresponding to the shaded histograms in 
figures [and [2] Now the aim is to find the six pairs of circles due to the dodecahedral 
topology. If there would be only the pure ordinary SW contribution, a value of Xmax = 1 
would clearly stand out of the distribution due to S?. However, the Doppler and ISW 
contributions degrade this signal to the grey histograms in figures [I] and [2] Thus 
the six pairs would probably yield values of Xmax between 0.4 and 0.8 for OQro¢ = 1.011 
and between 0.4 and 0.7 for Qot = 1.020, and values of U(p) between 0.3 and 0.75 
for both values of Qot, respectively. In the case of an analysis of a real sky map, one 
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would obtain many thousands of Xmax values which would not correspond to matched 
circles, i.e. having a probability distribution looking like that of S. Embedded in this 
distribution, one would for the dodecahedron obtain six Max values corresponding to 
real matched circles. Because of the large overlap of both distributions, the signal of the 
six pairs of circles could easily be swamped. Thus there is a non-vanishing probability 
that the six pairs of matched circles can be overlooked in the sky maps. The search for 
matched circles is in addition made more difficult by the finite thickness of the SLS and 
the Sunyaev-Zeldovich effect, a possible systematic in the removal of foregrounds and 
further secondary effects which are not included in our analysis. 


6. Discussion and Summary 


Taking the WMAP value 4 = 1.02 + 0.02 as its face value, we conclude that the 
Universe is slightly positively curved and therefore necessarily spatially finite, since all 
spherical space forms possess this property. At present, the only way to decide on the 
particular topology realized in our Universe consists in studying a concrete space form 
and then compare its predictions with the observational data. 

In this paper, we produce evidence in support of the hypothesis that the spatial 
structure of the Universe could be given by the Poincaré dodecahedral space D, provided 
that the total energy density parameter hot is in the range 1.016...1.020. For the 
detailed comparison with the WMAP data, it was crucial that we could base our 
predictions for the CMB anisotropy upon a large number of vibrational modes of 
the Poincaré dodecahedron. Since the eigenfunctions of the Poincaré dodecahedron 
are not known analytically, they have to be computed numerically. Using the very 
efficient algorithm described in section B] we could take in our calculations the first 
10521 eigenfunctions into account. In addition, we put forward the Conjecture (25) 
for the homogeneous spherical space forms. Combining the numerically computed 
eigenfunctions with the Conjecture, we were able to calculate the CMB multipoles in 
a broad range, see figures D] and [8] covering the region of the large scale fluctuations. 
Fixing the cosmological parameters as Qpar = 0.046 and OQmat = 0.28, h = 70, in 
agreement with the WMAP data, there are only two parameters left: the amplitude 
a of the scale-invariant Harrison-Zeld’ovich spectrum (88), and the density parameter 
Qrot > 1 or, equivalently, the curvature radius Ro = 7 (Oot — 1)-/?, (The dark energy 
is identified with a cosmological constant with magnitude Qa := Qiot — Qrad — Qmat-) 
The parameter aœ was obtained from a fit to the Cı values as determined by WMAP in 
the range l € [20,45]. (This differs from the normalization in [85], where a has been 
set to match the WMAP multipole with l = 4 exactly, since higher multipoles could 
not be calculated in due to the very limited number of available eigenfunctions.) 
Having fixed in this way all parameters except Qot, we studied the CMB multipoles and 
the temperature correlation function C(v) for the dodecahedron as a function of Qtot. 
One should keep in mind that we restricted the parameter space to the above values. 
Thus we cannot exclude that other choices could lead to different values for Qtot due to 
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geometrical degeneracies. 

As a quantitative measure of the quality of our theoretical predictions, we calculated 
the S statistic (8) for the angles p = 60° and p = 20°. The result is given as a function 
of Qtot in figure B] and shows for this statistic the remarkable feature of pronounced 
minima at almost the same Oto values, Qot € {1.016, 1.020]. In particular, figure B] 
shows that the Poincaré dodecahedral universe leads for hot in the above range to a 
suppression of the CMB anisotropy at large scales in agreement with the WMAP data. 
This is also demonstrated in figure), which shows the temperature correlation function 
C(¥), and in figures Bf) and Bh) showing the CMB power spectrum (Q44 = 1.019 is 
used in these figures). Figure Zr) demonstrates clearly that the dodecahedron model 
describes the WMAP data much better than the concordance model. At this point 
it is important to notice that the suppression of the CMB anisotropy at large scales 
occurring in the dodecahedron is neither obtained by fitting the theoretical curves to 
one of the first (suppressed) multipoles nor to the correlation function C(v). Rather the 
suppression is due to the finiteness of the Universe in combination with a small volume 
of the dodecahedron, V(D) = x = 0.16449..., which is nicely illustrated in figure 
where we compare the power spectra of the dodecahedron D and the simply connected 
S’, respectively. The structures seen in figures B]and [Bþ) in the CMB power spectrum 
at multipoles l < 20, are mainly caused by the ordinary Sachs-Wolfe (SW) contribution 
(dotted curve in Bb)) whose analytic expression is given in eq. (ZI). The observed 
structures are a direct consequence of the discrete mode spectrum of D possessing large 
gaps at small wave numbers 8, see eqs.(26) and (27%). 

The particularly striking suppression of the quadrupole and octopole is also 
exemplified by figure [7] which displays the probability distribution P(éT?) (l = 2,3) 
obtained from sky simulations. The probability distributions for the dodecahedron 
show pronounced peaks which overlap nicely with the lo band of the WMAP data. In 
contrast, the concordance model predicts much broader distributions producing large 
values for ôT? and dT}, seen also in figure E] 

Summarising the results discussed so far, it is fair to say that the Poincaré 
dodecahedral universe provides a natural explanation for the strange suppression of 
power at large scales as observed by COBE and WMAP, if hot is in the range 
1.016...1.020, which solves “the mystery of missing fluctuations”. At this point we 
would like to emphasise that the good agreement of the Poincaré dodecahedral universe 
with the WMAP data does by no means exclude that there are other manifolds which 
lead to similar conclusions. Rather, it may well be that there exists a “topological 
degeneracy” in addition to the well-known “geometrical degeneracy” (see e. g. [24] [25]). 
Since the value ot = 1.02 + 0.02 reported by the WMAP team depends on several 
priors and includes the lo deviation uncertainty only, there still exists the possibility 
of the Universe having negative curvature. In fact, e.g. the hyperbolic Picard space 
[26] 27] leads also to a suppression of the CMB fluctuations at large scales in agreement 
with the WMAP data. 

There remains to detect a more direct signature of the dodecahedral topology. 
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For this purpose, we have studied the circles-in-the-sky signature in section D] The 
dodecahedral universe possesses exactly 6 matched circles with radii 42.7° and 49.8° 
for Oro = 1.016 and Qot = 1.020, respectively. If there would be only the ordinary 
Sachs- Wolfe contribution, one would predict an almost perfect matching of the circles. 
However, as it is seen in figure B] there are important contributions from the integrated 
Sachs-Wolfe and the Doppler effect. The three contributions together destroy a perfect 
matching as is illustrated for a given circle pair in figure Dh). While the ordinary 
Sachs-Wolfe contribution displayed in figure Dþ) shows the expected matching, the 
Doppler and the integrated Sachs-Wolfe contributions shown in figures Dk) and DH), 
respectively, differ appreciably over the full angular range, and thus lead to the poor 
matching demonstrated in figure Dh). 

The quantity (p), eq. M9), has been introduced for a search of matched circles 
in the microwave sky. Figures [I] and [2] display the probability distributions of Xmax 
and (36°), respectively. One observes broad distributions for both the dodecahedron 
D and the covering space S?. While the distribution for D is peaked at higher values 
of Xmax than the distribution for S*, the two distributions have a large overlap. Due to 
the contamination by the Doppler and integrated Sachs-Wolfe effect, the distribution 
for D has no peak at Xmax ~ 1, but rather gives for the most probable Xmax values at 
which the 6 matching circles are expected to be observed the interval 0.4...0.8. 

In view of these complications, it seems that a search for matching circles in the 
first-year WMAP data is quite difficult, and thus the negative result reported in [8 
does not yet allow to conclude that the dodecahedral model is excluded. We hope that 
the second-year WMAP data and the CMB measurements from the Planck satellite 
will improve the situation considerably such that we shall be able to decide about the 
topology of the Universe in the near future. 

Let us close with a quotation from Einstein, who wrote in 1919 in describing the 
Universe of beings in a “finite and yet unbounded Universe” [61] [62]: “The great charm 
resulting from this consideration lies in recognition of the fact that the universe of these 
beings is finite and yet has no limits.” 
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